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. We discuss the properties of a one-dimensional lattice model of a driven system with two species 

» I ' of particles in which the mobility of one species depends on the density of the other. This model 

CIh, was introduced by Lahiri and Ramaswamy in the context of sedimenting colloidal crystals, and its 

■ continuum version was shown to exhibit an instability arising from linear gradient couplings. In this 
paper we review recent progress in understanding the full phase diagram of the model. There are 
three phases. In the first, the steady state can be determined exactly along a representative locus 
using the condition of detailed balance. The system shows phase separation of an exceptionally 
robust sort, termed strong phase separation, which survives at all temperature. The second phase 

■ arises in the threshold case where the first species evolves independently of the second, but the 
O ' fiuctuations of the first infiuence the evolution of the second, as in the passive scalar problem. The 

second species then shows phase separation of a delicate sort, in which long range order coexists 
with fluctuations which do not damp down in the large-size limit. This fluctuation-dominated phase 
ordering is associated with power law decays in cluster size distributions and a breakdown of the 
, Porod law. The third phase is one with a uniform overall density, and along a representative locus 

' the steady state is shown to have product measure form. Density fluctuations are transported by two 

kinematic waves, each involving both species and coupled at the nonlinear level. Their dissipation 
. properties are governed by the symmetries of these couplings, which depend on the overall densities. 

^ ' In the most interesting case, the dissipation of the two modes is characterized by different critical 

I exponents, despite the nonlinear coupling. 
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I. INTRODUCTION 



Much recent progress in the general area of nonequilibrium steady states has come about by focusing on simple 
driven lattice-gas models which have the same symmetries as the more complex physical situations they represent. 
, For instance, the asymmetric exclusion process — a simple model of hard-core particles performing biased random 
\^ ■ walks — is representative of a class of current-carrying states and growth processes involving a single scalar field . 
. Now, there are several examples of physical processes which involve two coupled scalar fields for such situations, 

the corresponding lattice model involves two sets of Ising variables, with the time evolution of each species depending 
on the instantaneous configuration of the other. Such systems can display interesting, new types of steady states, and 
L«J ■ in this paper we discuss their occurrence within a recently introduced lattice model of coupled fields. 

I This model was initially proposed by Lahiri and Ramaswamy (hence referred to as the LR model) in connection 
■ with collective effects in sedimentation, the process by which heavier particles settle in a lighter fluid. Interactions 
' between sedimenting particles are mediated by the fluid , and give rise to a coupled dynamics of two fields, namely 

, the concentration (the particle number density) and the tilt field (the orientation of the principal axes of the particle 
'"^ ' distribution). A distinctive aspect of the inter-field coupling in this case is that within a hydrodynamic description, 
^ ] the coupling involves linear terms, in addition to more customary nonlinear coupling terms. This essential aspect 
O is captured by the 1-d LR lattice model. Linear couplings can also arise in other contexts, for instance in reduced 
models of magnetohydrodynamics drifting flux lattices in superconductors p2[ , and the motion of stuck and 
flowing grains on a sandpile [|l^. The essential aspect of linear coupling is captured by the 1-d LR lattice model. 
In the sedimentation problem, analysis of the linearized hydrodynamic equations points to an instability in a certain 
^ , regime, but leaves open the question of the ultimate steady state reached by the system. The lattice model allows 
' this question to be answered; the instability can be studied by varying the transition rates within the model, and the 
resulting steady state can be analysed and characterized 0,0- 

The instability in question is towards large-scale clustering, and the question arises to what extent the phenomenon 
resembles phase separation, familiar from equilibrium statistical mechanics. In fact, the nature of phase separation 
in nonequilibrium steady states is an area of great interest at present ]l5[-jl7t , and the study of the LR model proves 
quite illuminating in this context. The unstable phase turns out to be strongly phase separated, much as in the 



three-component model of Evans et al |18 1^. The phase separation is robust, and survives at all temperatures. 
The behaviour is quite different, however, at the threshold of the instability, which corresponds to the particularly 
simple situation in which one of the two fields evolves autonomously while the other is driven by the fluctuations 
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of the first, as in the passive scalar problem. There is phase ordering in this case as well, but it is quite delicate. 
Fluctuations of the order parameter are very strong, but do not destroy the ordered phase rather, they change 
its character in an essential way, leading to the breakdown of familiar characteristics such as the Porod law for the 
decay of scaled correlations. Finally, moving away from the threshold, even the disordered regime of the model has 
some surprises to offer. In this regime, the density is homogeneous on a macroscopic scale, but fluctuations are 
transported ballistically by two sets of kinematic waves |^ . The waves are coupled nonlinearly, yet under special 
symmetry-determined conditions (i.e., on a special locus in parameter space) the dissipation of each is governed by a 
distinct critical exponent [p3| . This is the first such instance to come to light in a system with nonvanishing nonlinear 
coupling. In the kinematic wave phase in general, the equations are a one-dimensional reduced model for the dynamics 
of small fiuctuations of a fiux lattice drifting through a type II superconductor and also serve as a generalised 
Burgers equation for magnetohydrodynamics (MHD) pl| ]. 

In the remainder of this section, we discuss the occurrence of the instability within a linearized hydrodynamic 
description, and then define the LR lattice model and briefly discuss its phase diagram. In the sections that follow, 
we discuss respectively strong phase separation (in the unstable regime), fluctuation-dominated phase separation (at 
the threshold), and kinematic waves and their decay (in the stable regime). 



A. Hydrodynamic Description and Instabilities 

The hydrodynamic equations of a coupled-field system involve the density fiuctuations Spi^2ix, t) = ^1,2(2^, t) ~ Pi 2 
of the two species in question. When dealing with conserved fields, the equations governing their time evolutions 
may be derived from the continuity equations dpi/dt + dJi/dx = and dp2/dt + dJ2/dx — 0, where Ji{Spi,Sp2) 
and J2{Spi,Sp2) are the currents of particles of species 1 and 2 respectively. We may proceed by writing each of the 
currents Ji and J2 in terms of a systematic part 2*, a diffusive part —Ddpi^2{x,t)/dx, and noise 771^2- The coupling 
between the two species arises from the systematic parts, for each of J^^'' and J2^'^ depend both on 5pi and 5p2- We 
obtain the requisite equations on expanding in powers of the density fiuctuations. We find it convenient to write these 
equations in terms of the integrated density fields hi,2{x, t) = Sp{x' , t)dx' . Then, to second order, we obtain 
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Let us first examine the effect of keeping only the linear first derivative terms on the right hand side of the above 
equations. These terms, which are of primary interest to us here, were absent on symmetry grounds in coupled-field 
problems involving directed polymer motion and surface growth, considered respectively in references and [^,^. 

Let the eigenvectors of the 2x2 matrix I '^^^ '^^^ ) be e+ and e_, and let the corresponding eigenvalues be c+ and 



There are three cases to consider: 



C21 C22 



(A) Stable: If c+ and c_ are real (i.e. if A = (cn — 022)^ + 4ci2C2i > 0), they represent the speeds of two waves. 
The two waves involve the eigenvectors e+ and e_, each of which is composed of a linear combination of hi and 
/i2. To proceed, rewrite Eqs. |^ in terms of e+(x,t) and e~{x,t). Notice that in the rest frame of each wave, 
the other wave moves with a finite speed |c+ — c_|; it is impossible to make both stationary by a Galilean shift. 
The nonlinear terms describe the couplings between the two waves which govern their dissipation. 

(B) Threshold: The discriminant A vanishes if one of the off-diagonal terms, say C12, is zero. At the linear level, this 
would imply that the field hi evolves autonomously, but (since C21 ^ 0) the time evolution of the 'passive scalar' 
/i2 is infiuenced by hi. There is a tendency for particles of type 2 to be driven together. In the corresponding 
lattice model, this system shows phase separation accompanied, however, by large-scale fluctuations. 

(C) Unstable: There is an instability if c+ and c_ pick up an imaginary part (i.e. if A < 0). The solutions of 
the linear equations then have an unbounded exponential growth of fluctuations. This instability signals the 
advent of a new state which is qualitatively different from the statistically homogeneous state assumed at the 
outset. By studying the corresponding lattice model in the unstable regime, we find that the system undergoes 
macroscopic phase separation of a robust sort. 
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B. Lattice Model and Phase Diagram 



The Lahiri Ramaswamy model ||^ is the lattice counterpart of the continuum equations discussed above, and 
analysing its phases helps us understand the behaviour of the system in each of the cases (A), (B) and (C) above. 

The LR model is defined in terms of two sets of variables {<Ti} and {Tj_i} which reside on two interpenetrating 
sublattices; the former occupy the integer sites and the latter the half-integer bond locations of a one-dimensional 
lattice with L sites. Each ct^ and t^_i is an Ising variable taking on values ±1. They represent discrete versions of 
the density and tilt fields in the sedimentation problem: If Ui is 1, there is a particle at site i, and if ai = —1, 
there is no particle (— ). The variable rj_ i = 1 and —1, implies two values / and \ of the local tilt respectively. A 
typical configuration of the full system is thus: +\ — / — / + \ — / + / + / + \—- 

Both sets of variables are conserved, i.e. '}2,<Ji and "Y^t^^i and the associated densities p° — + cri)l'^L and 
P2 — ^(H-Tj_ 1 )/2i are constant. The linear terms in the hydrodynamic equations lead us to consider a r-dependent 
local field which guides the tr-current and vice versa. Thus, for instance, the Kawasaki exchange dynamics of the 
adjacent spins Ui and ai^i occurs at a rate which depends on t^+i. The moves and the corresponding rates are 
depicted below: 
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The rates may be written succinctly as 



W{(Ji ^ cr,+i;r^+i) = Tj+i(CTj -cr,+i) 

w/ A (Pi +P-2) , (pi -P2) , 

In a simple correspondence to a hydrodynamic model, discussed in section IV below, the coefficients in the continuum 
description (Eq. |l]) depend on the rates pi,p2,ri,r2, as also on the overall densities p° and p2 of the two species. 
This correspondence allows us to demarcate each of the regimes (A), (B) and (C) discussed above in terms of the 
parameters of the lattice model. 

Throughout the rest of the paper, we will assume that ri > r2. Thus particles tend to move preferentially downhill, 
and holes uphill. Holding the ratio ri/r2 fixed, it is then interesting to ask for the effect of varying pi/p2- The result 
is summarized in the phase diagram of Fig. 1, for the half-filled case '''i = S = 0- We discuss the characteristics 
of each phase in brief. 

Pi > P2 (the light grey region in Fig. 1): In the corresponding contimmm equation, the diagonal coefficients cn 
and C22 vanish, while ci2 and C21 both have the same sign. Thus A is negative, corresponding to the unstable case 
(C). Since pi > p2, a local peak (/\) tends to become a valley (\/) if a particle resides on it, while a valley tends to 
convert into a hill if a hole sits on it. These moves act in concert with the motion of particles down slopes (determined 
by Ti > r2), to promote segregation of both kinds of spins. In fact, using the condition of detailed balance, the 



steady state can be determined exactly along the the line ri/r2 = P1/P2 shown in Fig. 1 |14|. The state shows phase 
segregation of a particularly strong sort, as discussed in section 2. 

Pi = P2 (the thick central line in Fig. 1): This is the threshold case, with C12 = leading to vanishing A. It 
describes a semi-autonomous problem in which the r's evolve independently of the tr's; they undergo a symmetric 
exclusion process. The tr's on the other hand obey the dynamics of Eq. 0. The system shows a particularly delicate 



sort of phase separation, which is accompanied by large-scale fluctuations in the steady state l20 
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^ P2/P1 (1,1) 



Figure l:The phase diagram of the LR model at half-filHng (^(7i — J2'''i-i/2 = 0). The strongly phase separated 
phase (SPS) is separated from the kinematic wave phase (KW) by the threshold, semi-autonomous line pi = p2 
where there is fluctuation-dominated phase ordering (FDPO). The steady state can be found exactly along the line 
r\jri — pi/p2 (using the condition of detailed balance) and along rxjri = P2/P1 (using the condition of pairwise 
balance). 

Pi < P2 (the dark grey region in Fig. 1): In this case C12 and C21 have opposite signs, so that A is positive. The 
tendency of particles to slide downhill is countered by the tendency of particle-rich valleys to become hills. The 
resulting steady state has a uniform macroscopic density and slope, while density fluctuations are carried by two 
kinematic waves, each corresponding to a mode which involves both u's and r's. The steady state can be found 
exactly along the line ri/r2 = P2/P1 shown in Fig. 1, by using the condition of pairwise balance [p^ . Interestingly, 
the symmetries of the mode evolution equations can be tuned by varying the overall densities, and it is thereby 
possible to achieve a situation in which the dissipation of the two modes have different dynamical exponents in spite 
of being coupled at the nonlinear level ||2^ . 



II. STRONG PHASE SEPARATION 



In the unstable phase of the LR model {pi > p2\ ri > r2, the region marked SPS in Fig. 1), there is a bootstrap 
effect at work: the tendency of cr particles to fall downhill into valleys works in tandem with the tendency of particle- 
rich regions of the r landscape to become valleys themselves. This results in the formation of ever deeper valleys and 
large-scale clustering of particles in them, leading ultimately to a phase separated state. 

A similar sort of phase separation has been shown to occur in a three species ABC model in which a cyclic rule 
determines preferential rates of exchanges of adjacent particles [|l8|,^. For the case with an equal number of A, B 
and C particles, the steady state was found exactly, using the condition of detailed balance. 

In the LR model too, the exact steady state can be found through a detailed balance argument in the case 

f^i = '^Ti = 0, ri/r2 = Pi/p2 (Fig. 1). This happens because although the problem is defined in terms of transition 
rates, it is possible to construct a Hamiltonian with respect to which detailed balance is valid. This Hamiltonian has 
the form 

N 

H^eJ^hkMak (4) 

k=l 

where /ifc{r} is a height field (the discrete analog of h{x) in Eq. |l|) defined by 

k 

hk{T} ^^Tj^l/2- (5) 

The steady state is thus described by an equilibrium Boltzmann measure, and a configuration {tr, r} has weight 
exp(— 7i{(T, t}/T) where T is the temperature. This can be seen by noting that the energy change on interchanging 
neighbouring a spins is 
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(6) 



while the change of energy under exchange of r spins is given by a similar expression with cr's and r's interchanged. 
Thus the ratio of Boltzmann weights of the configurations after and before an interchange of spins, say (Ti <-> Ci+i, is 
exp{—2AEa/T). This equals the ratio of forward to backward transition rates if the ratio e/T is related to the rates 
pi and p2 by 

Pi 

The ratio r2/ri of forward and backward transition rates involving the interchange of adjacent r's is also related to 
the same weight factor. Thus, provided ri/r2 = P\lpi holds, the steady state weight of a configuration is proportional 
to exp{—7{)/T . Notice that H. involves long-ranged coupHngs between spins ai and Tj-1/2. 

The ground state exhibits complete phase separation, corresponding to all cr = 1 particles being at the bottom of 
the potential well formed by phase-separated r spins. Pictorially, the state is 

-\-\-\-\+\+\+\+/+/+/+/-/-/-/ 

There are two interfaces across which there is a change of a spin values, and two more involving the r's. The spacing 
between the four interfaces is L/A, and each is completely sharp in the limit T — > 0. The effect of nonzero temperature 
is to smear out each of the four interfaces over a finite length. The interfacial profile can be calculated but for 
our purpose here it suffices to note that the interface width is of the order of T, and has basically no effect on the 
bulk of the spins. Thus in an infinite sample, phase separation survives at all T. This unusual behaviour can be 
traced to the fact that the energy grows super-extensively (~ (length)^), and thus dominates over the entropy at all 
temperatures. This allows the strictures against phase separation in equilibrium one-dimensional systems pq | to be 
avoided. 

When the filling is changed away from ai ~ Ti ~ 0, the condition of detailed balance no longer holds, and 
there is no description in terms of an effective Hamiltonian. Nevertheless, kinetic arguments can be given to show 
that phase separation survives even when the filling is changed. The extreme robustness of the phenomenon leads us 
to call it Strong Phase Separation. 

The approach to the steady state is exceptionally slow. It occurs through a coarsening process which involves thermal 
activation over large barriers, so we expect a coarsening length growing logarithmically in time. Logarithmically slow 
coarsening was found in the ABC model as well . 

In concluding this section, we mention a related model which has quite different behaviour. Recall that in the 
LR model, the mobility of each spin species depends on the local density of the other species. What if the mobility 
depends not on density, but a higher derivative of the density? A case in which the coupling is to the second derivative 
of the density can be analysed through detailed balance considerations, but this time, the Hamiltonian has short- 
ranged interactions p^ . There is thus no phase ordering in this model at nonzero T, in contrast to the strong phase 
separation exhibited by the LR model in the unstable regime. 



III. FLUCTUATION-DOMINATED PHASE ORDERING 



In the marginal case pi = p2, ri > r2 (the line FDPO in Fig. 1), the evolution of the r's proceeds independently of 
the cr's. The fiuctuating r landscape provides a source of nonequilibrium noise on the it's which tend to fall downwards 
along local slopes and cluster in the valleys. Since the strong landscape fluctuations never cease, the clustering is 
never as complete as in the case of strong phase separation. Rather, as we will see below, the character of long-range 
order is strongly modified by fiuctuations, hence the appellation Fluctuation-dominated Phase Ordering. 

Starting from an initial random arrangement of the a particles, the extent of clustering can be estimated quite 
simply [ po| : In any fluctuating surface characterized by a dynamical exponent z, there are rearrangements of the 
profile over length scales ~ t^^^ in time t. This sets the scale for the base lengths of new valleys which form in time 
t. In the case at hand, where the scaling properties of the surface are those of the Edwards- Wilkinson model |^^, we 
have z = 2. This valley base length should set the scale for spatial clustering of particles, so we expect the equal time 
correlation function to follow the scaling form 

{ao{t)(Jo+r{t)) = C{r/C{t)) with Cr^t^/'. (8) 

This is confirmed by numerical studies. This sort of scaling behaviour is a general characteristic of phase ordering 
dynamics, which describes the evolution from an initially disordered state to one which is phase separated [p8|. 
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However, the scaling function here is quite different from that in usual cases. Normally, with a conserved order 
parameter, C{y) drops hnearly for small y = r/C{t), i.e. C{y) « Co(l — ci|?/|) as ?; — > 0. By contrast, we find C shows 
a cusp for small values of \y\, i.e. 

C(y)«Co(l-ci|y|") (9) 

with a ~ 0.5 | |20[ |. This is significant, since a linear dependence on |?/| (which goes under the name of the Porod law) 
is a simple consequence of having phases separated by sharp interfaces ||2^ . The breakdown of this law is the first 
pointer to the unusual character of the state that is developing. 

In steady state, C again shows a similar scaling form, with the coarsening scale C(t) being replaced by the system 
size L. The cusp for small argument remains unaltered. Turning to one-point correlation functions, a suitable measure 
|p9| for our conserved spin system is the magnitude of the Fourier components of the density profile 

Qik)^\^i2e^''n,\, fc = ^. (10) 

where Uj = (1 + (Tj)/2 and m runs over 1, L — 1. A signature of an ordered state is that in the thermodynamic limit, 
the average values {Q{k)) approach zero for all nonzero fc, but Q{k —>■ 0) remains finite. We numerically monitored 
these averages (• ■ •) over the ensemble of steady state configurations and found that {Q{k)) — » in the thermodynamic 
limit for fixed fc 7^ 0, but (Q*) approaches a nonzero value, where Q* = Q{k = ^). 

(Q*) provides a quantitative measure of phase separation. We find that (Q*) ~ 0.18 indicating that the steady 
state is ordered, though not perfectly (a perfectly ordered state has Q* ~ 0.32). As a function of time, Q* fluctuates 
strongly, reflected in the fact that the distribution P{Q*) remains broad, with a root-mean-squared deviation that 
remains finite even as L — » 00. One may ask whether phase ordering disappears when, in the course of its fiuctuations, 
the value of Q* falls close to zero. The answer is in the negative. We observe that a dip in Q* is accompanied by a 
simultaneous rise in the value of either Q{m = 2) or Q(m = 3). This implies that whenever the system loses a single 
large cluster (making Q* small) either two or three such clusters appear in its place (making the values of Q{m) for 
some small m go up). Thus the system remains far from the disordered state, and always has a few large particle 
clusters which are of macroscopic size ^ L. A numerical study shows that the average size of the largest particle 
cluster ~ 0.14L. 

We observe further that the particle and hole cluster size distributions in the steady state decay as a power-law: 

Prob{l) - (11) 

with ~ 1.8. This allows us to understand some features of the correlation functions discussed above. Within the 
independent interval approximation [^0| which ignores correlations in the lengths of successive clusters, the two-point 
correlation function can be calculated and shown to have the scaling form of Eq. || |^^. The cusp in the scaling 
function is reproduced as well, with the cusp exponent a being related to 9 hy a + 9 — 2 within this approximation. 
On a qualitative note, the power law cluster distribution makes the interfaces between phases broad and structured, 
so that the argument which leads to the Porod law does not hold any longer. 

The phenomenon of fiuctuation-dominated phase ordering, with the characteristics of broad distributions of order 
parameter fluctuations, power-law cluster distributions, and cusps in the scaled two-point correlation functions, can 
be shown analytically to hold for a related simple model of surface depth fluctuations. Both steady state properties 
and time-dependent coarsening correlation functions can be calculated within this model |2^j2^ , and the results bear 
out the picture of FDPO discussed above. 



IV. THE KINEMATIC WAVE PHASE 



In the stable phase of the LR model (ri > ^2 , p2 > Pi), particle motion and landscape fluctuations work oppositely, 
so that the resulting state is disordered. The density of both cr's and t's is uniform, but the dynamics of density 
fluctuations is interesting, and this is the aspect we focus on. 

To begin with, we note that the steady state can be shown exactly to have a product measure form, provided that 
ri = p2 and r2 = pi (the line marked Pairwise Balance in Fig. 1). To do this, choose new symbols to denote the 
values of a and r: Use 1 if a site or bond is occupied by a or / and for — or \. Then the elementary moves 
reduce to 100 001 and Oil 110 (with rate ri), and the reverse moves 001 100 and 110 — > Oil (with rate 
r2). Now, in steady state the total probability flux into each configuration C is equal to the total fiux out of C. A 
sufficient condition for this is the condition of pairwise balance [M: in-fluxes and out-fluxes must balance in pairs. 
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i.e. for every out-flux C C, there is a unique in-flux C" C such that W{C" C)Pss{C") = W{C C')Pss{C) 
where H^'s are transition rates, and Pss^s are steady state probabihties. Using intuition gained from particle hopping 
models, this condition can be shown to be valid |23) if Pss(C) — Pss{C ). Thus every allowed configuration is 
equally likely in steady state, which implies that P{C) is the product of independent site occupation probabilities 
in the thermodynamic limit. An immediate consequence is that multi-point correlation functions in the steady state 
decouple into products of one-point functions, allowing us to find the exact expression for the current of each species: 

Ji = (ri-r2K(l-p?)(l-2p°) (12) 

for the (T-particle current, and a similar expression with 1 and 2 interchanged for the r-particlc current. 

Hydrodynamic equations for density fluctuations can be derived using these expressions for the current, on following 
the procedure of Section 1.1. Further, by taking linear combinations of the fields, one can construct eigenmode fields 
h± which decouple at the linear level. These fields describe wave- like modes travelling with speeds c±. The time 
evolution of the field is governed by the equation 



dh+ dh, ^d'^h, ^ , . . . ^^,dh< 

— " ' D — h Quadratic Nonlmearities — K[- 



dt dx dx^ dx 



^ dx ' ^ dx ' 



-ri+{x,t) 

(13) 



and a similar equation holds for h-. The quadratic nonlinearities in this equation are of the same form as in Eq. |^. 
The coefficients of these nonlinear terms, as also the coefficient K of the cubic nonlinearity in Eq. |l^, depend on the 
densities p° and P2- The new noise terms ri± are also delta-correlated. While the fields h+ and /i_ are decoupled at 
the linear level, they are coupled through the nonlinear terms, so that each mode can infiuence the dissipation of the 
other. 

Some of the coefficients in Eq. vanish for certain choices of densities p° and Then, special symmetries arise 
in the evolution equations and the dynamical exponents associated with the wave modes may change. Four cases 
arise. 

(a) RI symmetry: Invariance under up-down refiection (R) h — > —h and under inversion (/) of space x — > —x. 

(b) RI symmetry: Invariance under h — > —h and not under x —>■ —x. 

(c) Rl symmetry: Invariance under x —^ —x and not under h —h. 

(d) RI symmetry: Invariance neither under x —x, nor under h — > ~h. 

A term like obeys RI symmetry. Terms like (^) and (§|)'^ obey RI symmetry. The RI symmetry is respected 
by the term (§j)^, while a term like (^) added to it breaks that and gives rise to RI symmetry. 

To illustrate the occurrence and effects of different types of symmetries, we consider three special pairs of densities 

(I) For Pi — P2 = ^: Eq. ^3 and its partner for /i_ reduce to a pair of coupled equations, with linear first and second 
derivative terms and cubic gradient nonlinearities. Quadratic nonlinearities are absent. 

The equations describe two kinematic waves moving with speed c+ = r /2 and c_ = — r /2 where r = ri — r2. The 
nonlinear couplings imply that each wave influences the evolution of the other. In order to study the dissipation of 
say the -f mode, it is essential to move to the frame which co-moves with it, through the Galilean shift a; — > x -f c^t, 
t t. 

Evidently, in this frame, the —mode has a speed c_ — c+. The evolution equation is invariant under —h^, 
h_ —h_ but not x — > — x, because of the linear and cubic nonlinear terms. The RI symmetry holds in the 
rest frame of mode. Similarly the dissipation of the —mode can be studied by going to a frame which co-moves 
with it. It is easily seen that RI symmetry holds in this frame as well. 

Each of the -I- and —modes thus have the same symmetry as a single scalar fleld evolving through Edwards- Wilkinson 
dynamics with additional cubic nonlinearities. For the corresponding single-fleld system, we have z — 2 (though with 
multiplicative powers of logarithms) ]3l[ | , and we might expect the same behaviour here. 

(II) For pi — 1/2 and P2 ^ 1/2, on going to either of the frames in which the -I- mode or the — mode are at rest, we 
see that the RI symmetry applies for each of the fields; the R symmetry is broken by by quadratic nonlinear terms, 
and I is broken because of linear first order and cubic terms. The most relevant terms at the linear fixed point are 
the quadratic nonlinear terms. Thus we would expect these terms to govern the dissipation and give rise to the KPZ 
value z — 3/2 |^ for both the modes. 

(III) For Pi — P2 ^ 1/2, an interesting situation arises. In the co-moving frame of the —mode, the pair of equations 
have symmetry under — > — /i_ and /i+ /i+ but not under x — > —x. Thus the h- field has Rl symmetry, while 
the moving h+ field has RI symmetry. The same symmetries hold in the rest frame of the -t- wave. Based on these 
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observations, we expect z = 2 for the —mode (perhaps with multiphcative logarithmic corrections), and z = 3/2 for 
the + mode. 

Thus based on simple considerations of symmetry, we would expect different types of behaviour in the cases I, II 
and III. These expectations arc borne out by both numerical simulations and analytic calculations. The equal-space 
unequal-time autocorrelation function F{t) ^ t^^ for the lattice model was investigated numerically for the three cases 
discussed above. Our Monte Carlo results for /3 are in accord with expectations based on symmetries of the contimmm 
equations. In case I, we find /3 = 1/4 with evidence for multiplicative logarithmic corrections for both modes. In case 
II, we find /3 = 1/3 for both modes, indicating KPZ behaviour. Finally, in case III, we find that one mode has type I 
behaviour, while the other one has type II behaviour, as expected on grounds of symmetry. 

An analytical treatment of the generalised Langevin equations (13) is straightforward to carry out |23). The 
renormalised propagators and correlation functions for and h_ are calculated in a 1-loop self-consistent scheme. 
In case III in particular the character of the perturbation theory is controlled by the presence of kinematic waves. 
In contrast to what happens in [ pl| , their effect cannot be shifted away by Galilean transformations, and we find as 
expected from the symmetry arguments above that z+ = 3/2 and z_ = 2. The calculations for the other cases (I and 
II) likewise give results consistent with the symmetry arguments above. 



V. CONCLUSION 



We conclude by pointing out some connections to other models and discussing some open questions. 

As discussed in Section 2, the strongly phase separated phase of the LR model is closely related to the corresponding 
phase of the permutation-symmetric ABC model ]l^ , pj| especially at the symmetric point of each model where the 
detailed balance condition is valid. However, an analog of the threshold case of fluctuation-dominated phase ordering 
which occurs in the LR model has not yet been found in the ABC or related models. In particular, the permutation 
non-symmetric three-component model of Arndt ct al p3t shows a transition from a phase-separated phase to a 
disordered phase with a large correlation length ||3^ ], but the state at the transition point does not seem to show 
FDPO. 

The threshold case of the LR model is tantamount to the problem of hard-core particles sliding down a fluctuating 
surface. This leads to the question of whether fluctuation-dominated phase ordering survives if the surface in question 
evolves differently. A study of the problem with KPZ dynamics for the surface shows that FDPO survives |2|], though 
the value of the cusp exponent changes from that in the EW case. On the other hand, with the surface evolving 
according to Das Sarma-Tamborenea dynamics [^5| , there is no cusp in C (y) (Eq. |^) , and the Porod law is recovered 
p2[ . Another interesting direction to investigate is what happens if the falling particles are noninteracting. Some 
results in this direction are available from a study of domain- wall dynamics in a model of growth of binary fllms . 

The continuum equations of motion in the wave phase are a one-dimensional reduced model of those governing the 
dynamics of small fluctuations of a flux lattice drifting through a clean type II superconductor [ p^ . These equations, 
or more properly the equations implied by (|^) for dxhi, dxh2, with some of the couplings set to zero, are Burgers-like 
equations for magnetohydrodynamics (MHD), as shown in [ pl| . The dynamics in this kinematic- wave region of the 
phase diagram is largely governed by KPZ-like critical exponents (in particular z = 3/2). On special loci in the 
space of parameters, however, weak dynamic scaling prevails: if h± are the eigenmodes of the linearised versions of 
the equations, then the complete nonlinear theory on these loci is governed by z+ = 3/2, z_ = 2 for and h- 
respectively p^ . 

Finally, an important open question concerns the behaviour of the LR lattice model in higher dimensions. Recently, 
it has been shown that a higher-dimensional version of the ABC model continues to show strong phase separation 
pTf . It would be interesting to investigate analogous questions for all phases of the LR model as well. So far, the 
only result along these lines is the demonstration of product measure along particular loci in the disordered phase in 
higher dimensions p^ . 
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